p-band in a rotating optical lattice 
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We investigate the effects of rotation on the excited bands of a tight-binding lattice, focusing 
particularly on the first excited (p-) band. Both the on-site energies and the hopping between 
lattice sites are modified by the effective magnetic field created by rotation, causing a non-trivial 
splitting and magnetic fine structure of the p-band. We show that Peierls substitution can be 
modified to describe p-band under rotation, and use this method to derive an effective Hamiltonian. 
We compare the spectrum of the effective Hamiltonian with a first principles calculation of the 
magnetic band structure and find excellent agreement, confirming the validity of our approach. We 
also discuss the on-site interaction terms for bosons and argue that many-particle phenomena in a 
rotating p-band can be investigated starting from this effective Hamiltonian. 



I. INTRODUCTION 

Ultra-cold atom experiments display amazing versatil- 
ity and promise to improve our understanding of many 
particle physics. There is now hope for direct experimen- 
tal realization of many models which were constructed as 
effective models of condensed matter systems, such as the 
Hubbard model [1]. However, the extent of the ultra-cold 
atom experiments are not limited to previously discussed 
models, novel systems such as dipolar 2] and spinor gases 
Q are created as well. Through the interplay of these 
experiments and theories aimed at explaining or stimu- 
lating them, a better understanding of quantum many 
particle physics emerges. 

One of the problems that has been discussed exten- 
sively in the condensed matter literature, but never ex- 
perimentally realized is the effect of a periodic potential 
under high magnetic fields When the magnetic flux 
through the unit cell of the periodic potential is of the 
order of one flux quantum, the energy spectrum displays 
complex magnetic fine structure within the Bloch bands. 
Conventional condensed matter systems with lattice con- 
stants in the order of nanometers require thousands of 
tesla magnetic fields to be able to see these effects, which 
is very far from being experimentally feasible. 

Recently a number of authors have argued that this 
model can be realized in a cold atom set-up within cur- 
rent experimental capability H H, 0, H, S OH [H[ • As the 
atoms are neutral the magnetic field is expected to be cre- 
ated either by rotating the optical lattice, or by optically 
induced potentials. A weak rotating optical lattice has 
already been realized in vortex pinning experiments [l2| , 
and there is ongoing work about creating effective mag- 
netic fields via light induced coupling [131 ]. Theoretical 
studies of these systems promise interesting phenomena 
such as lattice quantum Hall effects [3, [IH, [l|| [l?} or 
the observation of topological conductance quantization 
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Most of the recent theory, as well as the previous inves- 
tigations in the condensed matter literature are focused 
on the lowest (s-) band of the lattice. An important rea- 
son for this focus is that the magnetic fine structure of 
this lowest band is very well described by the Peierls sub- 
stitution [l9l |. The resulting band structure for the split 
s-band is easily obtained from a difference equation and 
is a self-similar fractal, known as the Hofstadter butterfly 
The original tight-binding Hamiltonian is modified 
only by the addition of phases to hopping amplitudes and 
serves as a starting point for the investigation of many- 
body physics in this system. 

An exciting development in cold atom physics has been 
the realization that higher bands in an optical lattice are 
also experimentally accessible [2(| Hl[ . The physics of the 
first excited band, the p-band, contains surprises such as 
Bose condensation at non-zero momentum [22jl . or or- 
bitally ordered Mott insulators [H [H, EE EH]. For a 
system of fermions p-band physics can be accessed triv- 
ially by filling the s-band completely; surprisingly the 
relaxation time for bosons in the p-band is long enough 
to allow experimental access to pure p-band physics. 

A natural question to ask about the p-band physics is 
how the particles in the p-band respond to the effective 
magnetic field created by rotation. Experimentally, if 
a strong rotating optical lattice is realized, the p-band 
should be as accessible as the s-band. One can ima gine 
the already rich physics of the p-band [IE EE l29l. l30l. 
l3ll . [32l [33l | to be strongly affected by the magnetic field, 
as both the orbital order within each lattice site and the 
hopping between different lattice sites will be modified. 
Beyond the single particle physics, it is not clear how the 
various many particle phases, such as orbitally ordered 
Mott insulators, will be affected by rotation. 

The theoretical investigation of such effects requires a 
consistent method of incorporating the phases generated 
by the magnetic field into the lattice Hamiltonian. For 
the s-band, Peierls substitution, in which one builds an 
effective Hamiltonian by replacing k with (p— e A /c)/h in 
the energy band function gives a satisfactory description 
of the one particle physics [4j . Starting from this effective 
Hamiltonian interaction effects can be investigated. The 
accuracy of Peierls substitution for the s-band has been 
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checked by numerical solutions of the Schrodinger equa- 
tion [19(. However, as for degenerate bands (of which the 
p-band is the simplest example) the conjecture was that 
"wherever the unperturbed Bloch bands touch or overlap, 
it is not possible to obtain the magnetic sub-structure by 
semiclassical methods, even approximately, by means of 
a universal rule for the whole Brillouin zone" [l9[. 

In this paper, we generalize the Peierls substitution 
procedure to the p-band, and obtain an effective Hamil- 
tonian for the p-band of the rotating optical lattice. We 
show that after an appropriate diagonalization in k space, 
which assumes temporarily that only the on-site ener- 
gies are affected by the degeneracy lifting field, Peierls 
substitution is still a good option to obtain the detailed 
magnetic fine structure. We check the spectrum obtained 
from the effective Hamiltonian with an accurate numeri- 
cal solution of the two-dimensional Schrodinger equation 
and obtain excellent agreement. This method should in 
principle be applicable to other degenerate bands and it 
provides us with a means to examine inter-particle inter- 
actions. 

The paper is organized as follows: In the next sec- 
tion, we introduce the Hamiltonian for a rotating optical 
lattice, and discuss the tight binding limit. Section III 
contains a discussion of the Peierls substitution scheme, 
the resulting magnetic fine structure, and its compari- 
son with direct numerical solutions. In Section IV, we 
give the effective Hamiltonian including interactions and 
conclude in section V. 



II. THE MODEL 



We assume that Vq is deep enough for a tight-binding de- 
scription to apply to the system and furthermore concen- 
trate on the dynamics of the particles in the first excited 
(p-) band of the lattice. Our approach is to first cast this 
Hamiltonian into a second quantized form which includes 
the anisotropic hopping between nearest neighbor sites, 
the on-site zero point energies, and also the shift caused 
by rotation. Not only do we expect the hopping between 
lattice sites to be affected, as it was for the s-band, but 
also the on-site energies to be modified. However, since 
the hopping and on-site Hamiltonians do not commute, 
a common transformation that accounts for both modi- 
fications cannot be found. 

To overcome this difficulty, we temporarily assume 
that the hopping amplitudes are not affected by the effec- 
tive magnetic field and the only change is in the on-site 
energies. Our expectation is that in this way we will ob- 
tain two non-degenerate bands to which we can apply 
Peierls substitution. This procedure is rather ad hoc the 
validity of which is later checked through a comparison 
with the first-principles results presented previously [l9[ 
and reproduced here partially. 

We proceed with considering the following p-band 
tight-binding Hamiltonian (the energy spectrum of which 
is measured relative to the center of the tight-binding s- 
band) for non-interacting particles including the on-site 
zero-point energies and the rotation term (— SIL Z ) [22l.l33j 

H = _E V(^ R+Q e, 5 ,R + H - C -) + ; -E^R 
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We start with the Hamiltonian for a particle in the 
rotating frame of a two-dimensional square lattice 
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+ Vb[sin 2 (fca;) 



+ -mLu\r 2 — fiz • r x pj_ 
iky)], 



(1) 



where = (p Xl p y ) and r = (x, y). m is the mass of 
the particle, u>± is the transverse harmonic trapping fre- 
quency, fl is the rotation frequency, and Vq is the depth 
of the optical potential created by a laser beam with wave 
number k = 2ir/\ (for counter-propagating laser beams 
lattice constant a is equal to A/2). In what follows, we 
use photon recoil energy E-r = h 2 k 2 /(2m) as the energy 
unit. This Hamiltonian can be rearranged as 



(pi - m!lz x r) 2 ,- 9/ 9/ 
H = — z + Vb[sin 2 (fca;) + sin 2 (fcy)] 
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+ ^m(Lu'i - n 2 )r 2 . 
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We neglect the last term assuming that fi is very close 
to lo±_ , so essentially we deal with a particle under an ef- 
fective magnetic field B = 2mc£l/e in a lattice potential. 



where the summation is over all lattice sites R and band 
indices y, = x,y (since the problem is two-dimensional, 
p z orbital will not be considered). As usual, b' R (& M p_) 

is the creation (annihilation) operator for a particle in 
the p M band at lattice site R, e v is the unit vector along 
the v direction, to is the frequency of the isotropic har- 
monic oscillator potential which models the lattice poten- 
tial around its minima, and i M „ is the anisotropic hopping 
amplitude. The explicit expression for t^ v (in the absence 
of rotation) is 
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-fr 2 V 2 



'"""" L 2m 
= tuS^u - (1 - d^tj 



^( r )]<?W( r + ae v )dr 



(4) 



where V(r) is the periodic lattice potential and 4> p (r) 
is the localized Wannier function corresponding to the 
Pfi band. When we approximate the lattice potential by 
a harmonic oscillator around a minimum, these can be 
expressed as a product of harmonic oscillator eigenfunc- 
tions, i.e. Pa .(r) = ui(x)u Q (y) and cf> Py {r) = u (x)ui(y), 

u n (x) being the rfi 1 harmonic oscillator eigenfunction. 
iy is the hopping amplitude between two neighboring 
p orbitals aligned along the orbital orientation and t± 
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FIG. 1: Lowest three bands for the 2-dimensional sinusoidal 
lattice potential. The energy difference between the lowest 
two bands (s and degenerate p levels) (measured from the 
band centers) is huu within the harmonic oscillator approxi- 
mation for the potential minima, uo being the oscillator fre- 
quency. For Vo = 2QE R , huj = 7.7739£ fl . 



is the amplitude when the orbitals are oriented trans- 
versely with respect to the line connecting them. Both 
amplitudes are defined to be positive and tii ^> t± due 
to larger overlap. Since the lattice potential is separa- 
ble in x and y coordinates, ty and t±_ indeed have simple 
expressions in reference to the one-dimensional problem. 
t± and iii are one quarter of the widths of the lowest 
and next lowest bands for V = Vo sin 2 (fcir), respectively. 
By solving the Schrodinger equation numerically, we find 
t ± = 0.0025£ fl and t\\ = 0.0603£ fl for V = 20E R . The 
on-site zero-point energy fno also has the simple inter- 
pretation of being the energy difference between s and 
p levels (bearing in mind the harmonic description, see 

Fig. mi. 



III. PEIERLS SUBSTITUTION AND 
MAGNETIC FINE STRUCTURE 

We perform a Fourier transformation on the Hamil- 
tonian [Eq. I[3"])] as a preliminary for diagonalization in 
momentum space. The transformed Hamiltonian is 



(5) 



where e ^ 



2 cos(fc„a). Since the Hamiltonian 

is bilinear in creation and annihilation operators, it is 
diagonalizable by a Bogoliubov transformation. Defining 
Ak = e xk. + ano - Ak = e yk + we observe that the 
Hamiltonian is diagonalized in k space by the following 



transformation: 
1 



a k = [( cos ^k + sin 6 k) b xk + 1 ( cos °k ~ sin 9 \d h v \s\ 

= 7^ [( COs6, k _ s[n6 k) b x\t _ *( C0S #k + S[n9 k) b yk\ >( 6 ) 



with 



cos 29^ = 



sin 2Q\r = 



'Ak_^k\ 2 ' 
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The diagonal Hamiltonian has the form 



^(k)^/3 k ], 



with 



Ak + Ak 



±nnyi + ( 



Ak - Ak\2 
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(J) 



where upper (lower) sign refers to a {(3). From this point 
on, we can apply Peierls substitution to the dispersion 
relation [Eq. I|7jl] to obtain an operator out of it, i.e. 
we change k to (p — eA/c)/ft using the Landau gauge 
A = Bxy. The resulting Hamiltonian is transparent 
only when expressed in terms of a power series 



Ak + Ak I toTi _i_ Vilt Ak i 
2 ±nu i i+ 2 [ 2HQ ' 

'Ak~ Ak\4 , i 



(8) 
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with the assumption that \ f lk ~ f 2k \/2Ml — \ cos(k x a) — 
cos(fc y a) | (t || +t±)/Ml is smaller than one. If (in +t±)/h£l 
is much smaller than one, terms of lower order in 
(/ik — Ak)/ 2 ^ wn ^ be more dominant and one needs 
to consider only few terms for a desired accuracy, in- 
stead of summing the whole series. Increasing accuracy 
is achieved by adding higher order terms. In a typical 
experimental condition, for instance, with Vq — 20Er 
and Ml ~ Er, the ratio {t\\ + t±)/Ml is ~ 0.063, so a 
first order approximation may be sufficient for the de- 
sired accuracy. Here, we give the results to second order 
m (Ak — Ak)/ 2 ^' ^ or completeness. The approximate 
energy band functions, where we retain terms up to sec- 
ond order, are then 



= + ci [cos(k x a) + cos(fcj,a)] 



: [cos 2 (fez a) + cos 2 (&J, a) — 2 cos(fc x a) cos(fc y a)] , (9) 



where 



h(u> ± Q), c\ = tit — t±, and cf = ±(t 



t±) /2m. After converting cosines into sums of exponen- 
tials and making the Peierls substitution we obtain dis- 
crete translation operators, which allow us to express the 
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eigenvalue problem as a difference equation. Since trans- 
lations along y are multiplied by phases depending on x 
in the Landau gauge, one should be careful in creating 
an operator from cross terms such as exp(ik x a) exp(ik y a). 
The correct way of transforming should yield Hermitian 
operators and is obtained by symmetric combinations 
such as 



a plane wave [H 

1>(x,v)=e***g{x). (11) 



„ip x a/h i(p y —eBx)a/h _i_ i(p y —eBx)a/h ip x a/h 
ik a ik a /1 r\\ 

e x e y — > . (1UJ Making the substitutions x = na and y = la, n and 

I being integers, and acting the effective Hamiltonian 
Due to the translational invariance of the problem along E at p[(p — eBxy/c)/h] on the wave function [Eq. (fTTj) ], 
y direction, the y dependent part of the wave function is we get the following difference equation 



[g(n + 2)+g(n-2)] + 



± 



Ci c; 

y ~ "2 

,± 



2 [cos(27m0 — k y a) + cos(27r(n + 1)0 — k y a)\ jg(n + 1) 

Ci c? 1 
— 1- [cos(27rn0 — k y a) + cos(27r(n — 1)0 — k v a)\ >g(n — 1) 



cos(47rn0 — 2k y a) + c\ cos(2-7m0 — fe^a) + c; 



T 1 '2 



$(n) = S fl (n), (12) 



where c^, ci, and were introduced following Eq. (|9j) 
and = a 2 B/(hc/e) is the magnetic flux quantum per 
unit cell. can be expressed in terms of the rotation 
frequency Q, as = 2ma 2 il/h. 

When = p/<?, p and q being relatively prime inte- 
gers, the difference equation [Eq. (fT2|) ] yields q equations 



together with the Bloch condition g{n + q) = e x9a g(rt) 
due to the q-site translational invariance in the x direc- 
tion. By diagonalizing the resulting q x q coefficient ma- 
trix for several k x and k y pairs, we obtain the energy 
eigenvalues which are plotted in Fig. [5] as a function of 
0. We observe that each split band further divides into 
q sub-bands forming a pattern which has close resem- 
blance to the Hofstadter butterfly. This result is in fact 
anticipated since = ±(tii + 1±) 2 /2hQ is much smaller 
than c\ = tii — t± and if we simply neglect it as a first 
approximation, the energy band function [Eq. ([9])] will 
just be that of the tight-binding s-band, except that we 
have Ml which gives rise to increasing separation between 
the split p-bands with increasing 0. Our approximation 
becomes poorer as (or f2) becomes smaller since we 
require that (t\\ + t±)/fi£l be small. This is apparent in 
Fig. [5] in which we highlight the region where two bands 
overlap. However, if we increase the lattice depth, which 
decreases the hopping amplitudes, we can increase the re- 
gion of validity. Equivalently, we can say that our results 
should improve as increases. Another improvement 
option would be to consider a higher order expansion in 
translation operators, which models long-range hopping 
with yet smaller amplitudes. 

To be able to judge the accuracy of the magnetic fine 
structure obtained by our method we compare it with 



FIG. 2: Magnetic fine structure of the p-band for Vb = 20Er. 
Two-fold degenerate zero-field p-band is split into two as = 
p/q grows. Each split band further has q sub-bands. Our 
approximation fails in the shaded region, corresponding to 
< 1/6 (the nearest = 1/q to 1/5, for which the spectrum 
is displayed in Fig. where two bands overlap. This region 
can be made narrower if the lattice depth Vo is increased. 



a direct numerical solution of the Schrodinger equation, 
starting from the Hamiltonian [Eq. ((2}]. One method of 
numerical solution is to reduce the problem to a magnetic 
unit cell using magnetic translation symmetry and solve 
the two dimensional Schrodinger equation within this 
unit cell using finite difference methods. Unfortunately, 
the magnetic unit cell size increases with q, the denomi- 
nator of the flux = p/q, and the non-trivial boundary 
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FIG. 3: (a) Approximate energy levels, corresponding to 
4> = 1/5, in our effective Hamiltonian approach, (b) Band dia- 
gram obtained through a first-principles calculation in which a 
truncated basis of harmonic oscillator wave functions is used. 
Dashed lines show the edges of the zero-field p-band. 



conditions required by magnetic translation symmetry 
makes this direct solution method computationally inef- 
ficient. Another, more efficient method, which was first 
developed byZak [34J, and then expanded on by Ober- 
mair et. al. [19( , is to use magnetic translation symmetry 
to reduce the two dimensional Schrodinger equation to a 
set of p one-dimensional equations with non-local cou- 
plings. This equation can be handled with relative ease 
using a truncated basis of harmonic oscillator wave func- 
tions. Still, a numerical calculation is efficient only for 
pure cases with <j) = 1/q and for small q values. 

In Figs. [3] and |U we compare our results with those 
obtained by a direct numerical calculation along the lines 
of Ref. [10| . Calculations with the effective Hamiltonian 
are much faster and the results are as good as the direct 
numerical solution. For instance, in the case of <fr = 1/5 
(Fig. [3]) the agreement is already good, but if we increase 
(j> to 1/3 (Fig. 2]), apart from a slight overall shift, we 
see that band gaps are also more faithfully reproduced. 
The computational efficiency of the effective Hamiltonian 
method for the single particle problem is striking, but its 
real utility is that it can be used as a starting point to 
include interactions in the system. 



IV. EFFECTIVE HAMILTONIAN 

Until now, we have essentially been dealing with the 
single particle spectrum. The results we obtained can be 
utilized to examine the case of many particles, if we first 
write the effective Hamiltonian in real space 



FIG. 4: Energy levels for (f> = 1/3. (a) Results of the ef- 
fective Hamiltonian approach, (b) First-principles band dia- 
gram. Our approximation is better compared to the case of 
4> — 1/5, depicted in Fig. [3j in the sense that here band gaps 
are also more correctly captured, apart from a slight over- 
all shift. Also shown, by dashed lines, are the edges of the 
zero- field p-band. 
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r and r' have x = na; 
otherwise. 



Here, (r, r') denotes summation over nearest neighbors in 
the square lattice (with separation a), ((r, r')) over next- 
nearest neighbors (with separation V2a), and (((r, r'))) 
over next-next-nearest neighbors (with separation 2a); 
± sign refers to the hopping direction. We note that the 
next-nearest and next-next-nearest coupling amplitudes 
turn out to be the same in our approximation. This ef- 
fective Hamiltonian represents non-interacting particles 
moving in the p-band of a square lattice under a particu- 
lar magnetic flux </>. The connection between the new and 
old operators is made through the following definition 



H, 



4 «<r,r'»> 



"4/3r<) 



cos + sin ^k&j. 



/3 k = cos6» k 6 k 



sin 6» k 6^, 



(14) 
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with 6^ = (b^ ± ibyfc )/V%- The operator (6, ) an- 
nihilates a particle with momentum Kk whose z compo- 
nent of angular momentum is —h (h). To first order in 
(i|| + t±)/Ml, and (3-^ are of the following form 



6f± 



t\\ +tj 
2hfi 



(cos k x a — cos k y a)b^, 



where the upper (lower) sign refers to a (/3). Af- 
ter expressing cosines as exponentials, we make the 
Peierls substitution, i.e. we change k to k — eBxy/hc 
in the coefficients of an d interpret the resulting 
factors exp(±i27t(/)x/a) as momentum translation op- 
erators whose action on a function of k is given by 
exp(±i2TT<fix / a) f (k) = /(k =F 2ir<frZ./a). Fourier transfor- 
mation of these modified operators yields the real space 
operators as 



*H + 1± 
AMI 

i27T$ni.=F 
e °n,Z+l 



+ 1,1 

— il-K&TL 



*L-i)> (is) 



where the indices (n, i) specify the a; (= na) and y (= la) 
coordinates. 

For bosons, the short-range repulsive interactions be- 
tween particles can be incorporated into our model as an 
on-site interaction energy which can be writtem up to 
terms renormalizing the chemical potential, as [22j | 



Hi, 



U 



E 



J2r 

3h 2 



(16) 



U = 9 \4> Px Jr)\*dr, 



where — Y^n fyLr^r ^ s the boson number operator, 

L^r = —ifo(bl;rbyr — b yr b x r) is the z component of the 
angular momentum of a boson at site r, and g > is 
the short-range repulsive interaction strength. The inter- 
action Hamiltonian can be written in a microscopically 



more revealing way using n v 
tion nr = + n Y and L z ? 
interaction becomes 



2U 



(b$ 
—h(riY 



In this nota- 
n Y ). So the 



H in t = — V [{n+f + {n Y ) 2 +An+n Y ]. (17) 
6 r 

By adding H mt [Eq. ([17])] to F e // [Eq. (fl3])]. we obtain 
the Hamiltonian for interacting bosons in the p-band of 
a rotating optical lattice. 



V. CONCLUSION 

We considered how the degenerate excited bands of a 
tight-binding optical lattice are affected by the effective 
magnetic field created by rotation. Specifically consid- 
ering the first excited (p-) band of a two dimensional 



lattice, we pointed out that the magnetic field causes 
not only the hopping between different lattice sites to 
be modified, but also changes the on-site energies. We 
showed that once the modification of the on-site energies 
are explicitly taken into account, the Peierls substitution 
scheme can be used to obtain an effective Hamiltonian 
and the energy spectrum of the system. 

The energy spectrum contains not only the splitting 
of the two bands under the effective magnetic field, but 
also the fine structure forming a pattern similar to the 
Hofstadter butterfly. We compare the energies obtained 
from the Peierls substitution procedure with a direct nu- 
merical solution of the Schrodinger equation, and observe 
that our procedure matches the numerical solution to a 
very good accuracy. 

The effective Hamiltonian is obtained by using a se- 
ries expansion in the ratio of the hopping parameter t\\ 
(3> tj_) to Ml, which is a small parameter for tight- 
binding lattices except in the limit of very slow rotation. 
We carry out this expansion to second order and the re- 
sulting effective Hamiltonian contains hopping between 
all lattice sites that can be connected by traversing two 
links [Eq. ([T5])]. 

While we performed a second order expansion, it is 
instructive to display the effective Hamiltonian to first 
order in (t^ + t±)/Ml in terms of the original operators 
b^r- 

Heff = ^r^E [Ct,t>Q4)%{1 ~ 25 (r) . (r0 .) +H.C] 
<r,r') 

+ t J _t ± ^ C Yr [{b^b+ + (b Y )%,] 
{T,r>) 

+ V" [h(w + n)njt" + H(w - f2)nf ] 
r 

= E W&U+aeAre££ +ae " ^ +H.C.) 
+ E 6 lr^r + iMl 5>LV - b\ Y b xV ), (18) 
: ±l2 ™*, (r) x = (r% = na; 



CV r' — 



1, 



(r) x ^ (r%. 



This Hamiltonian incorporates the first non- vanishing ef- 
fects of rotation and can be used as an effective Hamilto- 
nian if (tu + t±)/Ml is not large. Indeed a recent preprint 
which appeared while this paper was in preparation uses 
this form as a starting point [33|. However, to investigate 
corrections for slower rotation one has to go to higher or- 
ders as in Eq. (fT3"l) . In Eq. (fP8|) , we display the vector 
potential A explicitly to express the gauge invariance of 
the effective Hamiltonian. Our numerical work was car- 
ried out using a higher order approximation [Eq. (|13[) ] 
which is also gauge invariant. 

In conclusion, we showed how Peierls substitution can 
be used for degenerate bands and checked its accuracy 
with direct numerical solutions. By investigating how 
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operator transformations are modified through Peierls 
substitution [Eq. (|f 5[) ] we derived a first order effective 
Hamiltonian in real space [Eq. (fT8|) ]. 

Going to the next order, we obtain a more accurate, 
but more complicated effective Hamiltonian, which dis- 
plays how higher order hopping is modified by the effec- 
tive magnetic field. Finally we also give the expression 
for on-site interaction for bosons in terms of the angu- 
lar momentum 'up' and 'down' operators. We hope that 
our results stimulate further theoretical and experimen- 



tal investigations of the p-band physics under an effective 
magnetic field. 
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